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Abstract 

Many-particle many-hole states near the magic number 20 have been 
investigated in the relativistic mean field formalism using the fixed config- 
uration method. Neutron particle-hole states in neutron-rich nuclei with 
TV ~ 20, i.e. 30 ' 32 ' 34 Mg,and 28 ' 30 Ne are studied to find out the ground 
state configuration. The ground state of 32 Ne,Mg as well as possibly 
34 Mg comes out actually as particle-hole states. Proton-neutron excita- 
tion across the shell gap at 20 in nuclei with TV — Z, i.e. 32 S, 36 Ar, 40 Ca 
and 44 Ti is investigated as possible origins of superdeformed configura- 
tion. Observed superdeformed bands in 36 Ar and 40 Ca can be described 
as many-particle many-hole states. 

1 Introduction 

One of the most interesting topics in recent years in the study of low energy 
structure of light nuclei is the existence of highly deformed states in near- 
closed shell nuclei at relatively low excitation energy. Superdeformed (SD) 
bands have been observed in light TV = Z nuclei 36 Ar[Tl|5] and 4O Ca0El 
as well as in 38 Ar 5 with the band heads situated around 5 MeV excitation 
energy. Studies on 32 Mg and the neighbouring neutron rich nuclei with 
TV ~ 20 showed very large B(E2) values for the 2^~ — > 0^~ transition 
indicating a large quadrupole deformation [S] |7| |5] The existence of 
highly deformed solutions near the magic number 20 is now understood 
to be a result of many-particle many-hole excitation across the shell gap. 
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From the spherical shell model perspective, particles are promoted from 
the sd shell to the £7/2 orbital across the shell gap at N, Z — 20. When the 
prolate deformation is large, the oblate deformation driving last orbitals of 
the sd shell are higher in energy than the low Q orbitals coming from the 
f-r/2 orbital. Thus a state may be formed with a large prolate deformation 
at a relatively low excitation energy by filling up the single particle levels 
from below. 

In neutron rich nuclei around 32 Mg, a pair of neutrons may be ex- 
cited across the N — 20 shell gap resulting in the 2hu> configurations. 
In Z — 10 — 12 and N = 19 — 22 nuclei, recent shell model calcula- 
tions |10l 1111 1121 ll.'-i I predict these 2hu> states to he energetically below 
the normal Oftu states. In the Nilsson model picture, the shell gap van- 
ishes at large deformation for neutrons. This so called 'island of inversion' 
where the 2?uj states lie below the normal states have been extensively 
investigated both theoretically and experimentally. Different mean field 
calculations |14l 1151 1161 117| differ in their predictions for the ground state 
configurations of different nuclei in this region. There is considerable 
disagreement among the different calculations regarding the boundary 
of the island of inversion. Another area where the calculations do not 
agree is whether the neutron and proton distributions are identical in 
these nuclei. Pritychenko et al. have measured the transition probabili- 
ties B(E2; 0+ -> 2+) in the nuclei 26 < 28 Ne and 30 ' 32 ' 34 Mg to study the role 
of these intruder configurations in an effort to determine the boundary. 
They conclude that the ground state of the nuclei 26 Ne and 30 Mg are well 
understood in terms of normal Ohu configurations. On the other hand 
Chiste et al. have performed inelastic scattering experiment to measure 
the charge and mass deformation of 30,32 Mg. They conclude that both 
the isotopes show similar type of structure involving 2hu configurations 
and there is no significant difference between the proton and the neutron 
distributions. 

The SD band El in 40 Ca is based on the J* = 0+ state at 5.21 
MeV. In 40 Ca, low lying states with large collectivity have been explained 
in terms of neutron-proton excitation across the shell gap at N, Z = 20. 
Fortune et al. |18l I19| have suggested the dominant configurations of 
the J n = + states at MeV, 3.35 MeV and 5.21 MeV to be Op - 0ft, 
4p — 4h and 8p — 8h, respectively. A full shell model calculation involving 
both the sd shell and the fp shell is still beyond the present day capabil- 
ities. A recent truncated shell model calculation |20| has confirmed that 
the SD band in 40 Ca is indeed based on 8p — 8ft configuration. Cranked 
Hartree-Fock-BCS method has been applied by Bender et al. to study 
the ND and SD states in 32 S, 36 > 38 Ar and 40 Ca[£T]. In Ref. cranked 
relativistic mean field theory without pairing was applied to study the 
highly deformed states in 40 Ca. Zheng et aZ. |22ll2*5| have performed fixed 
configuration deformed Hartree-Fock calculation to study the deformed 
states in a number of nuclei including 40 Ca and 44 Ti. The origin of the 
SD band is predicted to be an 8p — 8h configuration in all the theoretical 
investigations. Configuration dependent cranked Nilsson-Strutinsky and 
truncated shell model calculations were performed in highly deformed 
states in 36 Ar. All calculations agree that in 36 Ar, these states are ex- 
pected to be based on 4p — 8h configuration coming from the excitation 
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of two protons and two neutrons across the shell gap. 

The aim of the present work is to study these large deformation multiparticle- 
hole states. We have applied the fixed configuration deformed mean field 
approach, a method applied earlier in this mass region by Zheng et al. 
|22II23| using a Skyrme Hartree Fock model. We propose to use the same 
method using deformed Relativistic Mean Field (RMF) approach. This 
approach has been very successful in describing various features of nu- 
clear structure. Ground state properties like binding energy, deformation, 
magnetic moment, isotopic shift and nuclear radius has been calculated 
with considerable accuracy [23] . It has also succeeded in describing excited 
states like normal deformed (ND) and SD bands, giant dipole resonances, 
etc. One of the prominent successes of RMF is the reproduction of the 
spin-orbit splitting naturally. It has been observed in |22| that the spin- 
orbit interaction plays a major role in the calculation for the excited states. 
The RMF approach naturally accommodates a correct spin-orbit interac- 
tion and thus, is expected to give a good description of the multiparticle- 
hole states. We have not come across any RMF calculation using the fixed 
configuration approach. 



Relativistic Mean Field theory is well known and has been described in 
detail elsewhere |24| . We have employed the force NL3 25 which was ob- 
tained by fitting the binding energy, charge radii and neutron radii of a 
number of spherical singly and doubly closed nuclei. This was explicitly 
designed for treatment of variation in isospin and has proved to be very 
useful in describing the ground state properties throughout the periodic 
table. Some of the nuclei studied in the present work, e.g. 32 ' 34 Mg and 
30 Ne are away from the stability valley and should be well described us- 
ing the force NL3. We assume axial symmetry and reflection symmetry 
and work using deformed harmonic oscillator basis. The method of cal- 
culation described in Refs. 20 ''27J has been followed. Throughout this 
work, 12 Fermion and Boson shells have been used for the calculation. 
The quadrupole charge deformation parameter j3c is obtained from the 
intrinsic charge quadrupole moment Qo using the prescription 



where _Ro = 1.2y4 1//3 , A being the mass number. The mass deformation 
parameter /3a can be similarly defined. 

The fixed configuration calculation is essentially very simple. The 
states studied in the present calculation are particle-hole states only in 
the sense of the spherical shell model. As already explained, at large 
deformation, the levels coming from the higher fp shell cross the upper 
levels of the sd shell giving rise to the configuration of these states. The 
level sequence obtained from Nilsson-like diagrams are used to construct 
solutions with large deformation. For example, in Fig. 1, the single par- 
ticle neutron orbitals in 40 Ca are plotted as a function of the quadrupole 



2 Calculation 




(i) 
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deformation parameter f3c- Single particle proton orbitals also show a 
similar kind of structure. The orbits are identified by the asymptotic 
quantum number Q n [NnsA] at large deformation and by Lj at zero de- 
formation. For N=20, there are two crossings. The first one is around 
(3 C = 0.3 between the 3/2+ [202] orbit and the 1/2" [330] orbit. The sec- 
ond crossing is at (3 C ~ 0.6 between the 1/2+ [200] orbit and 3/2" [321] 
orbit. In the fixed configuration approach, to obtain the 4p — \h solution, 
the 3/2 + [202] orbits are kept empty and the l/2~[330] orbits are kept 
filled for both types of nucleons. The relativistic Hartree equations are 
then solved iteratively to obtain the binding energy and the deformation 
self-consistently. The resulting solution is expected to have a deformation 
beyond 0.3. Similarly, the 8p — 8h solution obtained by further keeping 
the l/2 + [200] orbits empty and the 3/2~[321] orbits filled are expected to 
have a deformation beyond 0.6. The Hartree-BCS equations are then self- 
consistently solved to obtain the lowest energy solution. The solutions 
obtained in this way are of course intrinsic solutions, i.e. they contain 
contributions from all possible J values in the band. We are principally 
interested in the J n — + state, i.e. the band head of the K=0 bands. 
To obtain the energy of the + state we have used the cranking model 
formula 

E(0+) = E !n - (2) 

where the moment of inertia I cr is calculated using the cranking model 
formula |28| and j\ = J% + Jy. Here, Ei„ refers to the energy of the 
intrinsic solution. This is a very important correction as at a large defor- 
mation, it can be quite large while being zero for spherical configurations. 
For example, angular momentum projection in HF calculation in the su- 
perdeformed minimum of 40 Ca lowers the bandhead by 4 MeVl2T|. 

Pairing has been included in our calculation in the form of BCS method 
in the constant gap approach. We have taken both the proton and the 
neutron gaps to be equal to 1.0 MeV except for the closed shell configu- 
rations. In the light nuclei studied, the pairing energy is quite small. We 
have verified that our essential conclusions remain unaltered for reason- 
able variation of the gap parameters. 



3 Results 

We first present our results for neutron excitation across the N = 20 shell 
gap in a number of neutron rich nuclei near 32 Mg. Next we discuss the 
proton-neutron excitation across the shell gap N = Z — 20 in the N = Z 
nuclei, 32 S, 36 Ar, 40 Ca and 44 Ti. We present the results of our calculation 
in Table I. Here E(0 + ) refers to the energy of the + state belonging 
to the particular configuration and /3c and /3a refer to charge and mass 
quadrupole deformation parameters, respectively. 

The nucleus 32 Mg has N = 20 and Z = 12. In general, a ground state 
calculation without particle-hole excitation predicts the ground state to 
be spherical. This has been observed in earlier relativistic calculations 
also |16l I17| . In particular, although the deformation of a number of 
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light nuclei was calculated|17l I29| with some degree of accuracy in the 
relativistic hybrid derivative coupling model, the ground state of 32 Mg is 
still found to be spherical in this model. In the present calculation, we 
find that the spherical minimum has a binding energy 250.97 MeV. Next 
we study the 2p — 2h neutron excitation across the magic number 20. 
This corresponds to two holes in the z/3/2 + [202] state and two particles 
in the z/l/2~[330] state. The intrinsic solution has the binding energy 
248.93 MeV. After the correction due to angular momentum, the bind- 
ing energy of the deformed 2p2h + state comes out to be 252.50 MeV. 
Thus the v 1 — v~ 2 state becomes the actual ground state. The experi- 
mentally measured binding energy is 249.69 MeV |30| . The charge and 
the mass quadrupole deformation parameters come out to be 0.457 and 
0.492, respectively. The different experimental measurements on charge 
deformation parameters give the observed values /3c = 0.512 ± 0.044 [S] 
in intermediate energy Coulomb excitation and (3c = 0.61 ± 0.04 from 
inelastic scattering experiment [S|. From the results of another Coulomb 
excitation study T, a slightly lower value of (3c = 0.438 ± 0.046 can be 
derived. Thus the experimentally measured values are reasonably close to 
the theoretical result. The proton and the neutron distributions come out 
to be similar. This is in agreement with the measurements of Chiste et 
al. [Sj who have found no hint of any large decoupling between the proton 
and the neutron distribution. Our result also agrees with shell model cal- 
culations 1 1 II I13| as well as HFB results |15| . In contrast, another HFB 
calculation |I4| predict the ground state to be spherical, i.e. of normal 
configuration. We have also examined the u 4 — v~ 4 excitation but could 
not obtain any converged self-consistent solution. 

We have studied a few nearby even-even nuclei in the present formal- 
ism to find out the boundary of the island of inversion. In 30 Mg, the 
ground state turns out to be a Ohuj state with no p — h excitation. The 
binding energy of the ground state is calculated to be 243.66 MeV. The 
experimentally measured binding energy is 241.63 MeV 30 . The v 2 — v~ 4 
state is obtained by exciting two neutrons from the z^l/2 + [200] orbital to 
the !/l/2~[330] orbital across the N = 20 gap. The binding energy of the 
2p4h + state comes out to be 242.17 MeV, approximately 1.5 MeV above 
the ground state. The calculated charge and mass deformation values are 
0.522 and 0.601, respectively. Inelastic scattering gives the experimental 
value (3c = 0.52 ± 0.04 Coulomb excitation experiment gives a some- 
what lower value f3c = 0.43 ± 0.19|7|. As already mentioned, Pritychenko 
et aZ.0 have suggested that the ground state of 30 Mg to be a OTluj state 
while Chiste et aZ.[S| have concluded it to be a 2Tiui state. The energy 
systematics observed in our calculation seem to support the former alter- 
native although a configuration mixing calculation is required to reach a 
definite conclusion. 

In 34 Mg, the two solutions corresponding to the OTilo and 2hu> configu- 
rations come out to be very close in energy. The latter configurations has 
two holes in the !/3/2 + [202] orbital and two particles in the ^l/2~[330] 
and ^3/2~[321] orbitals each. The QTiuj state has binding energy 260.48 
MeV and the charge (mass) deformation values 0.345 (0.327). The bind- 
ing energy for the 2Tiu> state is 260.42 MeV and the f3c and (3a values are 
0.511 and 0.562, respectively. Coulomb excitation experiment provides 
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a upper bound for the charge deformation flc < 0.599 \7\. The experi- 
mental ground state binding energy is 256.59 MeV |30|. A more accurate 
treatment, particularly of the pairing correlation, may show that the 2?kj 
configuration is actually the ground state. 

In 30 Ne, the OTilo solution turns out to be spherical with the binding 
energy 215.93 MeV. The 2%uj 0+ state has the binding energy 216.62 MeV, 
making it the ground state. The experimental binding energy is 212.08 
MeV. In 28 Ne, the situation is similar to 30 Mg. The 2hu> state turns out 
to be 1.75 MeV above the OTkj state. The calculated binding energy of the 
ground state is 210.57 MeV while the measured value is 208.82 MeV. The 
neutron configurations of the 2fiu> states in 30 Ne and 28 Ne are identical 
with 32 Mg and 30 Mg, respectively. We have also studied 32 Ne but could 
not obtain any self-consistent solution corresponding to the v i ~v~ 2 state. 
We have also looked for low-lying 2Tuo states in 34 Si but could not find 
any self-consistent solution. 

Among the N = Z nuclei, we start with the doubly magic nucleus 
40 Ca. The binding energy of the ground state is calculated to be 341.92 
MeV, in good agreement with the experimental value 342.05 MeV. The 
ground state is spherical. A number of np — nh states has been studied 
in |22| under the SkHF formalism. They observed that for n=2 to 8, 
although (3 increases linearly with n, the solutions are nearly degenerate 
in energy. The 8p — 8h state has been identified with the SD band whose 
band head lies at 5.21 MeV excitation energy. The 4p — 4/i state has 
been identified with the ND band whose band head lies at 3.35 excitation 
energy. However for some of the Skyrme parameterizations, the + state 
of the ND band was found out to lie higher than that of the SD band in 
contradiction with experimental results. 

In the present calculation, we have concentrated on the 4p — 4/i ex- 
citation and the 8p — 8h excitation. The different configurations have 
been discussed earlier. In our calculation also, the ND band comes up 
very high in energy. In contrast, the excitation energy of the + state of 
the 8p — 8h band comes out to be 5.86 MeV. The transition quadrupole 
moment values also show a similar trend. The ND band quadrupole mo- 
ment is calculated to be 1.14 eb, large compared to the experimental value 
Qo — 0.74 ± 0.14eb[3"|. However, the SD band quadrupole moment comes 
out to be 1.89 eb, in excellent agreement with the experimental value 
Qo = 1.80lQ'2geb[j||. In contrast, the cranked RMF calculation without 
pairing 3, predicts Qo to be 2.0 eb. 

One of the problems in getting a good description of the ND band may 
be the near degeneracy of the np — nh states for even n. This degeneracy 
was observed in the present calculation also. A proper calculation should 
take into consideration the mixing between all the + states coming from 
the different np — nh configurations. This may provide a more accurate 
description. 

The ground state of the nucleus 36 Ar turns out actually to be oblate, 
This is in agreement with the number and angular momentum projected 
HF calculation of Bender et al. |21|. The calculated binding energy is 
306.36 MeV. In comparison, the experimental value is 306.72 MeV. For the 
4p — 8h state, we remove two protons and two neutrons from the 3/2 + [202] 
orbit and put them in the l/2~[330] orbit. The energy of the excited + 



6 



state after correction comes out to be 300.39 MeV. The excitation energy 
of the band head is thus 5.97 MeV. In comparison, this energy is calculated 
to be 5.90 MeV in Ref. |21| . The experimental value of the band head is at 
4.33 MeV. The experimental transition quadrupole moment is 1.18 ±0.09 
eb The calculated intrinsic quadrupole moment is 1.37 eb. 

We also study two other nuclei with N = Z, 32 S and 44 Ti where the 
SD band has not been observed. In 32 S, the experimental binding energy 
of the ground state is 271.78 MeV. We tried to look for solutions with 12 
holes in the 3/2+ [202], 1/2+ [200] and 5/2+ [202] orbitals and 4 particles 
in the 1 /2~ [330] orbitals. The excitation energy of the solution comes out 
to be 6.22 MeV. and the deformation to be very high, /3 ~ 1.0. We could 
not get any converged self consistent solution for only proton or neutron 
excitation. In 44 Ti, the experimental binding energy is 375.47 MeV while 
the calculated value comes out to be 372.48 MeV. We have obtained a 
8p-4/i solution with the 3/2+ [202] orbitals empty and the 1/2" [330] and 
the 3/2~[321] orbitals filled. The excitation energy of the corresponding 
band comes out to be 3.36 MeV and the charge deformation 0.530. In 
Ref. |22|. the excitation energy of the 8p — 4h band is found to be 5.6 
MeV and the deformation value (3 = 0.41. 

4 Conclusions 

We have observed that the nuclei 32 Mg, 30 Ne, and possibly 34 Mg belong to 
the island of inversion. As for the nuclei 30 Mg and 28 Ne, both nuclei with 
N = 18, the ground states are predicted to be moderately deformed Ohui 
states while the 2hu states have very large deformation. Measurements 
indicate that the ground state in these two nuclei are actually strongly 
deformed. The present calculation is unable to describe the ground state 
properties of these nuclei. In all the strongly deformed systems, we find 
that the proton and the neutron distributions are not decoupled, in agree- 
ment with Chiste et al. [Sj. 

In 40 Ca and 36 Ar, the SD band seems to be reasonably described 
in the present formalism. Particularly, the prediction of the intrinsic 
quadrupole moments is in excellent agreement with experimental mea- 
surements. However, in 40 Ca, the normal deformed configuration comes 
too high in energy. A source of error in the excitation energy of the 
bandhead may be the neglect of configuration mixing in the particle-hole 
bands. Surprisingly, we find no indication of any SD minimum in 32 S. 

We have included pairing through BCS approximation. The nuclei in 
the island of inversion have very large neutron excess and are close to the 
drip line. A more accurate treatment should involve the Bogoliubov ap- 
proach and possible coupling to the continuum as well as particle number 
and angular momentum projection. 

The calculations were done using the computer facilities provided un- 
der the DSA Programme by the University Grants Commission, New 
Delhi. 
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Table 1: Theoretical results obtained in the present calculation. See text for 
details. Configurations denoted normal are Ohco solutions. 



Nucleus 


Configuration 


E 


j'j in 2 


l rr /h 2 




0r 


0A 








MeV 




MeV -1 


MeV 








normal 




-250.97 






-250.97 


0.000 


0.000 




v 2 - v- 2 




-248.93 


33.67 


4.73 


-252.50 


0.457 


0.492 


30 Mg 

o 


normal 




-240.68 


10.95 


1.84 


-243.66 


0.285 


0.230 




v 2 - u- A 




-237.75 


36.93 


4.18 


-242.17 


0.522 


0.601 


34 Mg 

o 


normal 




-257.62 


21.53 


3.77 


-260.48 


0.345 


0.327 




i/ 4 - v- 2 




-255.81 


45.20 


4.91 


-260.42 


0.511 


0.562 


30 Ne 


normal 




-215.93 






-215.93 


0.000 


0.000 




v 2 - v- 2 




-213.06 


27.62 


3.88 


-216.62 


0.444 


0.502 


^«Ne 


normal 




-208.61 


5.32 


1.35 


-210.57 


0.202 


0.162 




v 2 - v-^ 




-204.72 


29.35 


3.58 


-208.82 


0.499 


0.619 


40 Ca 


normal 




-341.92 






-341.92 


0.000 


0.000 




V 2 Tt 2 - 1>~ 2 1T~ 


-2 


-328.43 


42.70 


5.50 


-332.31 


0.449 


0.442 




4 4 — 4 - 
V 7T — V 7T 


-4 


-329.80 


70.29 


5.61 


-336.06 


0.743 


0.733 


ae Ar 


normal 




-302.59 


12.20 


1.62 


-306.36 


-0.204 


-0.201 




is 2 ir 2 - v-^ 


-4 


-295.25 


47.65 


4.64 


-300.39 


0.642 


0.632 


:i2 S 


normal 




-266.12 


10.19 


1.43 


-269.68 


0.240 


0.236 




is 2 ir 2 - u~ 6 n 


-6 


-256.45 


64.52 


4.60 


-263.46 


1.030 


1.017 


44 Ti 


normal 




-370.80 


11.03 


3.28 


-372.48 


0.135 


0.133 




i/ 4 tt 4 - i/" V 


-2 


-364.45 


60.40 


6.46 


369.12 


0.530 


0.523 
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List of Figure captions 



Figure 1: The neutron single particle energy levels in 40 Ca as a function of 
quadrupole deformation j3c ■ The levels are indicated by the spherical quantum 
numbers at zero deformation and the asymptotic quantum numbers at large 
deformation. Continuous (dashed) lines represent positive (negative) parity 
levels. 
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